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EMCCD photometry reveals two new variable stars in the crowded 
central region of the globular cluster NGC 6981* 
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ABSTRACT 

Two previously unknown variable stars in the crowded central region of the globular cluster NGC 698 1 are presented. The observations 
were made using the Electron Multiplying CCD (EMCCD) camera at the Danish 1.54m Telescope at La Silla, Chile. The two variables 
were not previously detected by conventional CCD imaging because of their proximity to a bright star. This discovery demonstrates 
that EMCCDs are a powerful tool for performing high-precision time-series photometry in crowded fields and near bright stars, 
especially when combined with difference image analysis (DIA). 

Key words, globular clusters: individual: NGC 6981, stars: variables: general, stars: variables: RR Lyrae, instrumentation: high an- 
gular resolution 

1. Introduction 



* Based on data collected by MiNDSTEp with the Danish 1.54 m 

telescope A census of the va riable stars in the glo bular cluster NGC 698 1 

** email: skottfelt@astro.ku.dk was performed by iBramich et al.l (|201 lb (hereafter B 11). Using 

*** email: dbramich@eso.org data from 10 nights of observations with a conventional CCD 

^ email: rfiguera@eso.org; robertofiguera@gmail.com they were able to confirm the variability of 29 stars and refute 

* email: uffegj@nbi.dk the suspected variability of 20 others. Furthermore, 11 new RR 
^ Royal Society University Research Fellow Lyrae stars and 3 new SX Phoenicis stars were found in the 
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Study, bringing the total number of confirmed variable stars in 
NGC 6981 to 43. 

A problem with using a conventional CCD is that to obtain 
a reasonable signal-to-noise ratio {S /N) for the fainter objects 
in an image, the pixels in the brightest stars may well be satu- 
rated. For DIA, which is currently the best way to extract pre- 
cise photometry in crowded star fields (e.g. towards the Galactic 
bulge, in the central regions of globular clusters, etc.), the sat- 
uration of the brightest stars is even more problematic because 
the saturated pixels affect nearby pixels during the convolution 
of the reference image. Hence we cannot perform photometric 
measurements using DIA near saturated stars in conventional 
CCD images, which has a negative impact on the completeness 
of variability studies in crowded fields. In Bl 1, there are 4 sat- 
urated stars in the central region of their V reference image for 
NGC 6981, and it is therefore conceivable that their variable star 
census could have missed relatively bright variable stars (e.g. RR 
Lyraes). An improved completeness in variability studies makes 
it possi ble to draw fi rmer conclusions about Oosterhoff classifi- 
cation (ISmithllI995 ^ and to examine whether there is a gradient 
in the physical properties between the central and outer parts of 
the cluster 

Electron-Multiplying CCDs (EMCCDs) are conventional 
CCDs with an extended serial register where the signal is am- 
plified by impact ionisation before it is read out. This means 
that the readout noise is negligible when compared to the signal, 
even at very high readout speeds (10-100 frames/s), enabling 
the possibility of high frame-rate imaging. Numerous articles 
have described the possibility of using E MCCDs to obt a in ver y 
high spa tial resolution; see for instance Mac kav et al.l (l2004t) : 
iLaw et al. (2006). However, using EMCCDs to perform precise 
time-series photometry without throwing away photons (i.e. not 
Lucky Imaging) is a new area of investigation and the applica- 
tions are just starting to be explored. 

With high frame-rate imaging much brighter stars can be ob- 
served without saturating the CCD and the individual exposures 
can be combined into stacked images at a later stage in order to 
achieve the required S /N for the objects of interest. We note that 
EMCCD exposures need to be calibrated in a different way to 
conventional CCD imaging data. The algorithms required to do 
this ha ve already been developed and described bv lHarps0e et al.l 
(I2012h . 

Th e DIA technique, first introduced by lAlard & LuptonI 
( Il998l) . h as been i mproved by revisions to the algorithm pre- 
sented bv lBramich (2008); Bramich et"aD (l20 1 3 *) and is the opti- 
mal way to perform photometry with EMCCD data in crowded 
fields. This method uses a numerical kernel model instead of 
modelling the kernel as a combination of Gaussian basis func- 
tions and can thus give better photome tric precision even in very 
crowded regions (lAlbrow et al.ll2009i) . The method is also es- 
pecially adept at modelling images with PSFs that are not well 
approximated by a Gaussian. 

Using the superior resolution provided by high frame-rate 
imaging EMCCDs in tandem with DIA we can probe the sur- 
roundings of bright stars for variable stars which are inaccessi- 
ble with conventional CCD imaging. Using this technique we 
are able to present EMCCD photometry of two new RR Lyrae 
stars in the central region of NGC 698 1 . 

2. Data and Reductions 

The data were obtained over two half nights (26th and 27th 
August 2012) at the Danish 1.54m Telescope at La Silla 
Observatory, Chile, using the Andor Technology iXon-H model 



Table 1: Details of the two new variable stars found in 
NGC 6981. The celestial coordinates correspond to the epoch 
of the reference image, which is the heliocentric Julian date 
~2456167 d. The epoch of maximum light is given as a helio- 
centric Julian date (2450000 +) in column 5 and the period is 
given in column 6. 



Variable 


Var. 


RA 


Dec. 


^ max 


P 


Star ID 


Type 


(J2000.0) 


(J2000.0) 


(d) 


(d) 


V57 


RRl 


20 53 27.12 


-12 32 13.9 


6166.779 


0.334 


V58 


?" 


20 53 27.38 


-12 32 13.3 


6166.76 


0.285 



'"' We are unable to classify this variable (see Sec.|3) 



897 EMCCD camera. The imaging area of the camera is 
512x512 16;um pixels with a pixel scale of 0'.'09 which gives a 
45 x45 arcsec^ field-of-view (FOV). With such a small FOV, we 
chose to target the crowded central region of NGC 6981 includ- 
ing the saturated stars from Bl 1. The camera is equipped with a 
special long-pass filter with a cut-on wavelength of 650nm. The 
cut-off wavelength is determined by the sensitivity of the camera 
which drops to zero 0% at 1050 nm over about 250 nm. The fil- 
ter thus corresponds roughly to a combination of the SDSS i' -Hz' 
filters (.BesseIL2005) . A total of 44 observations with a frame- 
rate of 10 Hz were obtained. Each observation contains between 
3000 to 3500 exposures. 

Using the algorithms described in lHarps0e et al.l (l20I2l) . 
each exposure is bias, flat, and tip-tilt corrected, and the instanta- 
neous image quality (PSF width) is found. Then, for each obser- 
vation, the exposures are combined into images in two distinct 
ways: 

Quality-binned: exposures are grouped according to a binning 
in image quality and combined to produce images that cover 
a range in point-spread-function (PSF) width. 

Time-binned: exposures are grouped into time bins of width 2 
minutes to achieve a reasonable S /N at the brightness of the 
RR Lyrae stars. As opposed to Lucky Imaging, all frames are 
used. This gives a total of 125 data points in each light curve. 

To extract the photome try from the time-binned images w e 
used the DanDIA pipelin^] (lBramichl2008HBramich et al.i2013h . 
The pipeline has been modified to stack the sharpest of the 
quality-binned images to create a high-resolution reference im- 
age from which the reference fluxes and positions of the stars are 
measured. The reference image, convolved with the kernel solu- 
tion, is subtracted from each of the time-binned images to create 
difference images, and, in each difference image, the differential 
flux for each star is measured by scaling the PSF at the position 
of the star (see Bl 1 for details). Note that we have further modi- 
fied the DanDIA softwar e to employ the appropriate noise model 
for EMCCD data (Harps0e et al.ll20I2h . 

3. Results 

In order to detect new variable stars in our data, we constructed 
and visually inspected an image representing the sum of the 
absolute-valued difference images with pixel values in units of 
sigma. We found two new variable stars which we assign names 
V57 and V58, and the details of which are given in TablefT] Both 
stars are located close to a bright star as can be seen in Fig. [1] In 



' DanDIA is built from the DanlDL libraiy of IDL routines available 
at http://www.daiiidLco.uki 
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Fig. 1: Top: A cut-out of the B 11 V reference image correspond- 
ing to the field of view of the EMCCD camera. Bottom: Finding 
chart (constructed using our reference image) for the variables 
confirmed by Bl 1 and the two new variables V57 and V58 (note 
that V58 is not the star that is located at the upper right edge 
of the circle). North is up and East is to the right. The image 
size is 45 x 45 arcsec^. Notice the greatly improved resolution 
compared to the B 11 finding chart. 



the B 11 data, both of these variables are within the area that can- 
not be measured because of the saturated pixels from the bright 
star Using the saturation limits from Bll it can be concluded 
that the bright star is brighter than 14th magnitude in V. In our 
data we find that the bright star is about 4 magnitudes brighter 
than the RR Lyrae stars, which suggest that it is V~ 13 mag. 
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Fig. 2: Light curves, plotted in differential flux units, for the two 
new variable stars. Left and right panel show the first and sec- 
ond nights, respectively. The typical photometric uncertainty is 
plotted as an error bar in each panel. 
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Fig. 3: Phased light curves for the two new variable stars, plot- 
ted in differential flux units and using the periods from Table [1] 
Black and red dots represent the data from the first and second 
nights of observation, respectively. The typical uncertainty in the 
period is about 0.01 d for both variables 



The light curves for the two variables are shown in Fig. |2] 
There is increased scatter towards the end of each night which is 
due to a combination of high airmass and deteriorating seeing. 
The variable star periods were estimated using the string-length 
statistic Sq dP woretskvl 1 1 9 8 3h and the phased light curves are 
shown in Fig. [3] 

V57: With a period of 0.334 days, a sinusoidal-like light 
curve, and a brightness on the reference image similar to that 
of the other RR Lyrae stars, we can safely classify this variable 
as a first-overtone RR Lyrae star (RRl). 

V58: This object has no detectable PSF-like peak in the ref- 
erence image even though it shows clear PSF-like variations in 
the difference images. Hence the associated object is fainter than 
the cluster RR Lyrae stars. The period of 0.285 days is typical of 
an RRl star, although the light curve clearly deviates somewhat 
from being sinusoidal with relatively flat peaks and sharp drops 
in intensity. This star could be an eclipsing binary or an RRl 
star behind the cluster However, due to the lack of decisive evi- 
dence for either classification, we prefer to leave the variable as 
unclassified. 

The discovery of a new RRl variable in NGC 6981 changes 
the mean period of the RRl stars from 0.308 d (Bl 1) to 0.312 d. 
The updated ratio of the number of RRl to RR Lyrae stars is 
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found to be ~0.17 (compared to ~0.14 in Bll). Both of these 
quantities still agree very well with the classification of NGC 
6981 as an Oosterhoff" type I cluster (iSmithil995,) . 

4. Conclusions 

Using EMCCD data with DIA we found two previously un- 
known variable stars in the crowded central region of the globu- 
lar cluster NGC 698 1 . We have classified one variable as a first- 
overtone RR Lyrae and we have been unable to classify the other. 
The discovery of the new RRl star consolidates the classification 
of NGC 698 1 as an Oosterhoff" type I cluster 

Both variables are located in a crowded field and close to a 
much brighter star The previous study by Bll employing con- 
ventional CCD data with DIA failed to find these variables. Our 
discovery of these new variables in a carefully studied globular 
cluster is thus one of the first results to demonstrate the power of 
EMCCDs for high-precision time-series photometry in crowded 
fields. This means that EMCCDs can improve the results in a 
number of areas in astrophysical research, for instance the search 
for Earth-mass exoplanets in gravitational microlensing events, 
or, as mentioned here, a better contraint on the physical parame- 
ters of globular clusters. 
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